Optical conductivity of disordered graphene beyond the Dirac cone approximation 
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In this paper we systemically study the optical conductivity and density of states of disorded 
graphene beyond the Dirac cone approximation. The optical conductivity of graphene is computed 
by using the Kubo formula, within the framework of a full vr-band tight-binding model. Different 
types of non-correlated and correlated disorders are considered, such as random or Gaussian poten- 
tials, random or Gaussian nearest-neighbor hopping parameters, randomly distributed vacancies or 
their clusters, and random adsorbed hydrogen atoms or their clusters. For a large enough concentra- 
tion of resonant impurities, a new peak in the optical conductivity is found, associated to transitions 
between the midgap states and the Van Hove singularities of the main 7r-band. We further discuss 
the effect of doping on the spectrum, and find that small amounts of resonant impurities are enough 
to obtain a background contribution to the conductivity in the infra-red part of the spectrum, in 
agreement with recent experiments. 

PACS numbers: 72.80.Vp, 73.22.Pr, 78.67.Wj 



I. INTRODUCTION 



An important part of our knowledge on the elec- 
tronic properties of graphene, which consist of a two- 
dimensional (2D) lattice of carbon atoms,'i' can be de- 
duced from optical spectroscopy measurements (for re- 
cent reviews see Refs. '2'and'3). Infrared spectroscopy ex- 
periments allows for the contr ol o f interband excitations 
by means of electrical gatingjI^El similarly as electrical 
transport in field effect transistors. Within the simplest 
Dirac cone approximation, only vertical in wave-vector 
space transitions across the Dirac point are optically ac- 
tive, leading to a constant value for the optical conduc- 
tivity of undoped graphene of ctq = ire^ /2h. This leads 
to a frequency-independent absorption of -na ~ 2.3%, 
where a = /he « 1/137 is the fine structure constant. 
This fact was observed for suspended graphene in ex- 
periments in the visible range of the spectrum'^ and it 
was late r co nfirmed by further experiments in suspended 
graphen^^ and epitaxial graphene on SiC substrate.'^l 
For doped graphene with nonzero chemical potential /z, 
at zero temperature, in the absence of disorder and with- 
out considering many body effects, the allowed excita- 
tions are only those between particle-hole pairs with an 
energy difference larger than 2/j,, due to Pauli's exclusion 
principle. This would lead to a zero infrared conductivity 
below the energy w = 2/x, and the optical conductivity 
should be simply a step function a {uj) — ctoO (w — 2/x). 
However, a background contribution to the optical con- 
ductivity between < w < 2/i was observed in Refs. 5 
and [71 pointing out the relevance of disorder and many 
body effects. Another characteristic of the optical spec- 
trum is the Drude peak, which is built from a transfer 
of spectral weight from the low-energy interband conduc- 
tance to the w — )■ region of the spectrum,^^ although 



a strong suppression of the Drude peak at infrared ener- 
gies has recently been observed.l^ Furthermore, the flat- 
tening of the TT-bands at energies away from the Dirac 
point is responsible for the strong peak in the spectrum 
at higher energies (of the order of 5eV) which is asso- 
ciated to optical transit ions between states of the Van 
Hove singularities .'S'i^'i^l Finally, a method to control the 
intermediate excited states in inelastic light scattering 
experiments has been also reported, revealing the impor- 
tant role of quantum interference in Raman scattering.tSl 
This intense experimental work has been accompanied 
by a series of theoretical studies which have treated the 
problem of the optical conductivity at different levels of 
approximation.'iSHlll Por example, it has been suggested 
that the presence of spectral weight in the forbidden re- 
gion of the optical spectrum of doped graphene (below 
w = 2/i) can be associated to disorderj ^--^"^ electron- 
electron interactiorPSl or excitonic effects. I^Sl In particular, 
the effect of electron interaction in the spectrum has been 
considered in Refs. I27H35I Furthermore, understanding 
the role played by the different kinds of disorder that 
can be present in this material is essential to increase the 
mobility o f the sa mples. Besides the long-range charged 
impurities j ^^ l ^^ l '^ '' I other possible scattering sources such 
as ripples,'2S' strong random on-site p oten tials,!^ large 
concentration of hydrogen adatoms,'22rstrairP2El] or ran- 
dom deformations of the honeycomb lattice have been 
considered.'*^ 

In this paper, we perform a systemic study of the op- 
tical spectrum of graphene with different kinds of disor- 
der for both doped and undoped graphene, such as the 
randomness of the on-site potentials and fluctuation of 
the nearest-neighbor hopping. Special attention is paid 
to the presence of resonant impurities, e.g., vacancies 
and hydrogen adatoms, which have been proposed as the 
main factor limiting the carrier mobility in graphencl^^lf^ 
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Furthermore, depending on the way how the defects are 
distributed over the lattice sites, each kind of disorder 
can be either non-correlated or correlated. The non- 
correlated one corresponds to the case with uniformly 
random distributed disorder sources, i.e., the potential or 
hopping are randomly changed within a certain range, or 
the resonant impurities (vacancies or hydrogen adatoms) 
are randomly positioned over the whole lattice; the cor- 
related one means that the distribution of the disorder 
follow particular topological structures, such as Gaus- 
sian potentials or Gaussian hopping parameters, resonant 
clusters with groups of vacancies or hydrogen adatoms. 
In the present paper, we consider a noninteracting tt- 
band tight-binding model on a honeycomb lattice and 
solve its time dependent Schodinger equation (TDSE) to 
calculate the density of states (DOS). From this, the op- 
tical conductivity is calculated numerically by means of 
the Kubo formula. 

The paper is organized as follows. In Se c. [iTl we give 
details about the method. In Sees. Ull] and IIVI we show 
results for the optical conductivity of undoped graphene 
in the presence of non-correlated and correlated disorder, 
respectively. In Sec. |V]we calculate the optical spectrum 
of doped graphene. Our main conclusions are summa- 
rized in Sec. |Vll 



II. MODEL AND METHOD 

The tight-binding Hamiltonian of a disordered single- 
layer graphene is given by 

= - ^ {Uja\bj + h.c) + ^ Vic\ci, +H„np, (1) 

<i,j> i 

where where a] (bi) creates (annihilates) an electron on 
sublattice A (B), tij is the nearest neighbor hopping pa- 
rameter, Vi is the on-site potential, and Himp describes 
the hydrogen-like resonant impurities: 

H,mp^edJ24d^ + vJ2{dh^+^■c) , (2) 

i i 

where Sd is the on-site potential on the "hydrogen" impu- 
rity (to be specific, we will use this terminology although 
it can be more complicated chemical species, such as var- 
ious organic groupd^ and V is the hopping between car- 
bon and hydrogen atoms. For discussions of the last term 
see, e.g. Refs. 1391 Hll and |47l The spin degree of free- 
dom contributes only through a degeneracy factor and 
is omitted for simplicity in Eq. Q. A vacancy can be 
regarded as an atom (lattice point) with and on-site en- 
ergy — >■ 00 or with its hopping parameters to other 
sites being zero. In the numerical simulation, the sim- 
plest way to implement a vacancy is to remove the atom 
at the vacancy site. 

The numerical calculations of the optical conductivity 
and DOS are performed based on the numerical solution 



of the TDSE for the non-interacting particles. In gen- 
eral, the real part of the optical conductivity contains two 
parts, the Drude weight D {uj ~ 0) and the regular part 
(a; ^ 0). We omit the calculation of the Drude weight, 
and focus on the regula r part . For non-interacting elec- 
trons, the regular part ipSEH 

fafl (w) = lim — — / e~^*sinwt 

e^0+ Ojil Jq 

x2Im ((^1/ {H) J„ (t) [1 - / {H)] M^) dt, 

(3) 

(we put h = 1) where (3 = l/ksT is the in- 
verse temperature, is the sample area, / {H) ~ 
1/ l^e^'^'^^^-' + l] is the Fermi-Dirac distribution opera- 
tor, Jq [t) ~ e*^* JqC"'''^* is the time-dependent current 
operator in the a {— x or y) direction, and \lp) is a ran- 
do m sup erposition of all the basis states in the real space, 
i.e.JMHll 

W)=Y,a4\0), (4) 

i 

where are random complex numbers normalized as 
X^i kil^ ^ 1- The time evolution operator e~*'^* and the 
Fermi-Dirac distribution operator / [H) can be obtained 
by the standard Chebyshev polynomial representation.!^ 
The density of states is calculated by the Fouri er tra ns- 
form of the time-dependent correlation function^^MIl 

1 1"°° 

p{e)^- j J'Hv\e-'"'W)dt, (5) 

with the same initial state \ip) defined in Eq. Q. For a 
more detailed description and discussion of our numeri- 
cal method we refer to Ref. [511 In this paper, we fix the 
temperature to T = 300K. We use periodic boundary 
conditions in the calculations for both the optical con- 
ductivity and the density of states, and the size of the 
system is 8192 x 8192 or 4096 x 4096. 

III. NON-CORRELATED DISORDER 

A. Random on-Site Potentials or Nearest-Neighbor 
Hopping Parameters 

We first consider two different kinds of disorder: ran- 
dom local change of on-site potentials and random renor- 
malization of the hopping, which correspond to the diag- 
onal and off-digonal disorders in the single-layer Hamil- 
tonian Eq. Q, respectively. The former acts as a local 
shift of the chemical potential of the Dirac fermions, i.e., 
shifts locally the Dirac point, and the latter arises from 
the changes of distance or angles between the Pz orbitals. 
In order to introduce the non-correlated disorders in the 
on-site potentials, we consider that the on-site potential 
Vi is random and uniformly distributed (independently of 
each site i) between the values —Vr and -\-Vr- Similarly, 
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FIG. 1. (Color online) Numerical results for the density of states (left panels) and optical conductivity (right panels) of undopped 
graphene with different kinds of non-correlated disorders: (a,b) random on-site potentials, (c,d) random hopping parameters, 
(e,f) random distributed vacancies, and (g,h) random distributed hydrogen adatoms. Size of the samples is 4096 x 4096 for 
DOS, and 8192 x 8192 for optical conductivity. In the right column, the insets show a zoom of the optical conductivity in the 
infrared region of the spectrum. 
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the non-correlated disorder in the nearest-neighbor hop- 
ping is introduced by letting tij be random and uniformly 
distributed (independently of couple of neighboring sites 
between t — and t + tr- The presence of each 
type of disorder has quite similar effect to the density of 
states [see the numerical results with different magnitude 
of disorders in Fig. [T] (a) and (c) for the random on-site 
potentials {vr/t = 0.2, 0.5 and 1) and random hoppings 
[tr/t = 0.1, 0.3 and 0.5) respectively]. The spectrum 
is smeared starting from the Van Hove singularities at 
\E\ = t, and the smeared region expands around their 
vicinal areas as the strength of the disorder is increased, 
whereas the spectrum around the vicinal region of the 
neutrality point keeps unaffected unless the disorder is 
too strong. As the optical conductivity is proportional 
to the density of states of the occupied and unoccupied 
states, one expects a peak in the spectrum of the opti- 
cal conductivity at the energy cj sa 2t, which corresponds 
to particle-hole excitations between states of the valence 
band with energy E k, —t and states of the conduction 
band with energy E « tW These processes contribute 
to the optical conductivity with a strong spectral weight 
due to the enhanced density of states at the Van Hove 
singularities of the 7r-bands. Because we are considering 
a full TT-band tight-binding model for our calculations, 
this peak is also present in our results for the optical 
conductivity, as it is evident in Figs, [ijb) and (d) at 
Lo/t « 2, in qualitative agreement with recent experi- 
mental results. 1^ Notice that the height of the peak is 
sensitive to the presence of disorder, getting more and 
more smeared as the strength of disorder is increased. 
On the other hand, for this kind of disorder, for which 
there is no big change in the DOS around the Dirac point, 
one expects that the low energy spectrum of the optical 
conductivity should be robust for small disorder, i.e., the 
optical conductivity should follow the same spectrum as 
the clean sample without any disorder. These expec- 
tations are exactly what we observed in the numerical 
results of a (a;) shown in the insets of Fig. [l] (b) and 
(d). This is indeed the part of the spectrum that can 
be accounted for within the continuum (Dirac cone) ap- 
proximation. We can conclude that the non-correlated 
random disorder in the on-site potentials or hopping in- 
tegrals have almost no effect on the electronic properties 
(density of states and AC conductivity) in the low energy 
part of the spectrum unless the disorder is too large. On 
the other hand, the high energy inter-band processes be- 
tween states belonging to the Van Hove singularities of 
the valence and conduction bands are quite sensitive to 
the strength of these two kinds of disorder. 



B. Random Distributed Vacancies or Hydrogen 
Impurities 

Next, we consider the influence of two other types of 
defects on graphene, namely, vacancies and hydrogen im- 
purities. Introducing vacancies in a graphene sheet will 



create a zero energy mode (midgap state), effect that 
has been anticipated in many theoretical works, '^^ '^^^^ '^ 
and which has been recently observed experimentally 
by means of scanning tunneling spectroscopy (STM) 
measurements.'^ It is shown that the number of midgap 
states increases with the concentration of the vacancies^Sl, 
and the inclusion of vacancies brings an increase of spec- 
tral weight to the surrounding of the D irac poin t {E = 0) 
and smears the van Hove singularities! ^^ * '^^ ^ This is in 
fact the behavior found in Fig. [I] (e) for the DOS of 
graphene with different concentrations of vacancies rij,, 
where the numerical results with rix = 1%, 5%, 10% are 
represented and compared to the density of states of clean 
graphene. 

The presence of hydrogen impurities, which are in- 
troduced by the formation of a chemical bond between 
a carbon atom from the graphene sheet and a car- 
bon/oxygen/hydrogen atom from an adsorbed organic 
molecule (CH3, C2H5, CH2OH, as well as H and OH 
groups) have quite similar effect to the ele ctron ic struc- 
ture and transport properties of graphene .l^^l^ll xhe ad- 
sorbates are described by the Hamiltonian -ff^mp in 
Eq. ([1]). The band parameters V ~ 2t and td ~ — 1/16 
are obtained from the ah initio density functional theory 
(DFT) calculations.'*'^ Following Refs. [39] and EH we call 
these impurities as adsorbates hydrogen atoms but ac- 
tually, the parameters for organicgroups are almost the 
same.ESl As we can see from Fig. [l] (g) , small concentra- 
tions of hydrogen impurities have similar effects as the 
same concentration of vacancies to the density of states 
of graphene. Hydrogen adatoms also lead to zero modes 
and the quasilocalization of the low-energy eigenstates, 
as well as to a smearing of the Van Hove singularities. 
The shift of the central peak of the density of states with 
respect to the Dirac point in the case of hydrogen impu- 
rities is due to the nonzero (negative) on-site potentials 

Ed- 

The similarity in the density of states leads to simi- 
lar optical spectra for graphene with random vacancies 
or hydrogen adatoms, as it can be seen in Fig. [l] (f) 
and (h). In the high and intermediate energy part of 
the spectrum it is noticeable, apart from the smearing of 
the w « 2t peak due to the renormalization of the Van 
Hove singularities, the appearance of a new peak at an 
energy uj ^ t. This peak is associated to optical tran- 
sitions between the newly formed midgap states (with 
energy E ^ 0) and the states of the Van Hove singular- 
ities (with energy E ^ t). Notice that, contrary to the 
uj Ri 2t peak, the height of this uj t peak grows with the 
strength of disorder, due to the enhancement of the DOS 
at the Dirac point. Therefore, we expect that this peak 
should be observed in optical spectroscopy measurements 
of graphene samples with sufficient amount of resonant 
scatterers. 

In the low energy part of the spectra, the new struc- 
ture of the DOS around the Dirac point leads to a mod- 
ulation of the infrared conductivity, as it can be seen 
in the insets of Figs. [T] (f) and (h). The lower peaks. 
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FIG. 2. (Color online) Numerical results for the DOS (left panels) and optical conductivity (right panels) of undopped graphene 
with different kinds of correlated disorders: (a,b) Gaussian potentials, (c,d) Gaussian hoppings, (e,f) vacancy clusters, and (g,h) 
hydrogen clusters. The distribution of the clusters of impurities used for the results (e)-(h) are sketched in Fig. [s] 
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which in Figs. [T] (f) and (h) corresponds to a conductiv- 
ity a ~ 0.9(To for different concentration of impurities, 
might have their origin from excitations involving states 
surrounding the zero modes (central high peak in the 
density of states). At slightly higher energies there is 
a new set of peaks that can be associated to processes 
involving states at the boundaries of the midgap states. 
The optical conductivities in the region between these 
two peaks are in general smaller compared to those in 
clean graphene, which can be due to the fact that the 
midgap states are quasilocalized states. 

IV. CORRELATED DISORDERS 

A. Gaussian Potentials and Gaussian Hoppings 

As discussed in the previous section, the change of on- 
site potential can be regarded as a local chemical poten- 
tial shift for the Dirac fermions. If the random potentials 
are too large, characteristics of the graphene band struc- 
ture such as the Dirac points or the Van Hove singular- 
ities can disappear completely, and the whole spectrum 
becomes relatively flat over the whole energy range.l^ 
Therefore in order to introduce large values of random 
potentials but keep a relatively similar spectrum, in this 
section we use small con centra tions of correlated Gaus- 
sian potentials, defined a^^^ 

-^=Ec^^-p(^-fc|^), (6) 

where N^^^ is the number of the Gaussian centers, which 
are chosen to be randomly distributed over the car- 
bon atoms (rfe), Uk is uniformly random in the range 
[— A„, Ai,] and d is interpreted as the effective potential 
radius. The typical values of d used in our model are 
d = 0.65a and 5a for short- and long-range Gaussian 
potential, respectively. Here a « 1.42A is the carbon- 
carbon distance in the single-layer graphene. The value of 
TVj'J^p is characterized by the ratio P„ = N^^^/N, where 
N is the total number of carbon atoms of the sample. 
As one can see from Fig. [2ja), in the presence of locally 
strong disorders (A„ = 3t and t for short and long range 
Gaussian potentials, respectively) the whole spectrum of 
DOS is quite similar to the case of clean graphene, but 
with the emergence of states in the vicinal area around 
the Dirac point, and also a smearing of the Van Hove 
singularities. This kind of disorder leads to regions of 
the graphene membrane where the Dirac point is locally 
shifted to the electron {Uk < 0) or to the hole {Uk > 0) 
side with the same probability, rising the DOS at zero 
energy. The final spectrum is similar to the one of clean 
graphene but with a series of electron-hole puddles which 
are formed at the maxima and minima of the potential. 
The enhancement of the DOS around the Dirac point 
leads to the possibility for new excitations in the low en- 
ergy part spectrum, as compared to the clean case, as it 



can see in Fig. [2|^b). For the cases we consider, the pres- 
ence of long-range Gaussian potentials change the low 
energy optical spectrum completely with the emergence 
of a new peak around cj « 0.15i. The optical conductivity 
in the region uj < 0.24i is larger than in clean graphene 
but becomes smaller for ut > 0.24t. The increase of the 
conductivity might have its origin in the possible exci- 
tations between electron and hole puddles. Indeed, the 
renormalization of the spectrum obtained by consider- 
ing long-range Gaussian potentials leads to a larger opti- 
cal contribution than for short-range Gaussian potentials, 
which yield infra-red spectra much more close to that of 
a clean graphene membrane. 

The local strong disorder in the hopping between car- 
bon atoms is introduced in a similar way as the corre- 
lated potentials, i.e., with a distribution of the nearest- 
neighbor hopping parameter given bj^^ 

- i + E exp (-^^^^^^^) , (7) 

where iV*^^ is the number of the Gaussian centers (rfc), 
Tk is uniformly random in the range [— At,At] and dt 
is interpreted as the effective screening length. Simi- 
larly, the typical values of dt are the same as for the 
Gaussian potential, i.e., dt — 0.65a and 5a for short- 
and long-range Gaussian random hopping, respectively, 
and the values of iV*^^ are characterized by the ratio 

Pt = NLp/N. 

Numerical results for the DOS and optical conductiv- 
ity of graphene with short- (Aj = 3i, dt — 0.65a) and 
long-range (Aj = It, dt = 5a) Gaussian hoppings are 
shown in Fig. [2jc-d). This kind of disorder accounts 
for the effect of substitutional impurities like B or N in- 
stead of C, or local distortions of the membrane. Con- 
cerning the physics around the neutrality point, in this 
case the Dirac point remains unchanged although there 
is a local renormalization of the slope of the band. As a 
consequence, the Fermi velocity around Dirac is locally 
increased (when tk > 0) or decreased (when tk < 0). 
However, no midgap states are created by this kind of 
disorder, and the DOS remains quite similar to the one 
of a clean graphene layer, as it can be seen in Fig. ^c). 
In particular, the absence of an impurity band at w 
makes that the optical conductivity presents only slight 
deviations as compared to the clean case. This can be 
seen in Fig. [2]jd), where (apart from the smearing of the 
Van Hove peak) the optical spectrum, especially in the 
infra-red region, remains practically the same as in the 
absence of disorder. 



B. Vacancy Clusters and Hydrogen Clusters 

Correlated resonant impurities are introduced by the 
formation of groups of vacancies or adsorbed hydrogen 
atoms (see Fig. l3|. The center of the formed vacancy 



7 




FIG. 3. (Color online) Sketch of a gaphene sheet with vacan- 
cies (left pannels) or hydrogen adatoms (right pannels). The 
vacancies are presented as missing carbon atoms, whereas the 
hydrogen adatoms are highlighted in red. From top to bot- 
tom: Resonant impurites are distributed as the formation I 
(i? = 0), II (0 < 7? < 3a) and III {R = 3a) as described in the 
text. For illustrative purposes, the size of the sample shown 
in this sketch is 60 x 40, and the concentration of impurities 
is approximately equal to 2%. 



or hydrogen cluster (tc) is randomly distributed over the 
honeycomb lattice sites, with equal probability on both 
sublattices A and B. Each site (i) whose distance to one of 
the centers {R= \r — rd) is smaller than a certain value 
(i?c), is assumed to be part of the cluster, i.e., being a 
vacancy or adsorbing a hydrogen atom. We further in- 
troduce another freedom of the resonant clusters namely 
that their radius can change within the sample, allowing 
for a graphene layer with cluster of impurities of different 
size. This means that the value of Rc for each resonant 
cluster can either be different and randomly distributed 
to a maximum value, or can be kept fixed for all the clus- 



ters in the sample. We want to emphasize that as the 
center of the cluster is located on a particular sublattice 
A or B, the formation of the cluster does not preserve 
the sublattice symmetry and therefore can lead to the 
appearance of midgap states. 

Firts, in Fig. [2je) and (g), we compare the density of 
states with the same total number of resonant impurities 
(vacancies or hydrogen adatoms) but with different kinds 
of formations. We consider three different situations, i.e., 
randomly distributed uncorrelated single impurities (for- 
mation I), randomly distributed correlated clusters with 
varied radius of clusters (formation II) or with fixed ra- 
dius of clusters (formation III). The different structures 
are sketched in Fig. [3] Notice that the formation I is 
a limiting case of the formation III with all the radius 
of clusters being zero. As we can see from the results 
of the simulations, the number of midgap states is larger 
in the cases of uncorrelated single resonant impurities, 
and smaller for the case of fixed radius of resonant clus- 
ters. This is expected since the midgap states are states 
which are quasilocalized around the vacancies or carbon 
atoms which adsorb hydrogen atoms.'^^l^IIIIIIIl Therefore, 
for the same concentration of impurities, the number of 
midgap states will grow with the isolation of the impu- 
rities in small clusters. Something similar happens for 
the case of hydrogen clusters. This can understood by 
looking at Fig. |4j where we present contour plots of the 
amplitudes of quasieigenstates at the Dirac point or out- 
side the midgap region. The quasieigenstate |5'(e)) is 
a superposition of the degenerated eigenstates with the 
same eigenenergy e, obtained by the Fourier transforma- 
tion of the wave functions at different timeJ^ 

1 

'^te-*|^(i)), (8) 

where \ip {t)) = e^*^* \(p) is the time evolution of the ini- 
tial state \ip) defined in Eq. Although the quasieigen- 
state is not exactly the energy eigenstate unless the cor- 
responding eigenstate is not degenerated at energy e, 
we can still use the distribution of the amplitude in 
the real space to verify the quasilocalization of the zero 
modes in the presence of random impurities, "^^ or ob- 
tain the DC con ductivity at certain energies or carrier 
densities.'S^'^^'^ As we can see from Fig. [ij the con- 
tour plots of the quasieigenstates of graphene with va- 
cancy and hydrogen clusters are quite similar, i.e., the 
amplitudes on the carbon atoms which adsorb an hy- 
drogen atom are almost zero, just like if they are va- 
cancies. Furthermore, at the Dirac point (left panels of 
Fig. [ij corresponding to £' = 0) the quasieigenstates 
are semi-localized around the edge of the clusters (see 
the red color in the regions around the cluster). On the 
other hand, for energies above the impurity band, the 
states are not localized around the resonant cluster, and 
the amplitudes of the quasieigenstates are more or less 
uniformly distributed over the sample except within the 
clusters, where the amplitudes are zero. Therefore, as 
we have discussed above, for a given concentration of im- 
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FIG. 4. (Color online) Contour plot of the amplitudes of quasieigenstates at energy E = or E = 0.1094t. The radius of the 
resonant clusters is fixed at Rc = 5a. 



purities, the number of carbon atoms which are located 
around an impurity will be larger in the formation I than 
in the formation III. Then, the number of zero modes 
is also larger in I than in III, leading to spectra for the 
DOS and optical conductivity similar to the ones of clean 
graphene for samples in which disorder is concentrated 
in a small number of big clusters (formation III) than 
spread into a large number of small clusters (formation 
I), as it can be seen in Figs. [2]je)-(h). Finally, notice that 
the possibility for new excitations between the impurity 
and the carrier bands, leads to a modulation of the op- 
tical conductivity (as compared to the clean membrane) 
whose peak structure depends on the renormalized DOS 
and band dispersion of each case. 



V. OPTICAL CONDUCTIVITY OF DOPED 
GRAPHENE 

So far we have discussed the effects of disorder on the 
optical response of undoped graphene. In this section, 
we study the optical conductivity of graphene for finite 
values of the chemical potential, taking into account the 
effect of disorder. At zero temperature, a clean sheet of 
gated (doped) graphene has a zero optical conductivity 



in the region lu < 2/i, and an universal conductivity of 
<j{uj) = (To, due to optically active inter-band excitations 
throug h the Dirac p oint, for energies above the threshold 
^ > 2^j io | i5 | 2i | 29 l39l ^j^g presence of the disorder, the 

broadening of the bands as well as the appearance of pos- 
sible midgap states leads to a more complicated selection 
rule for the optical transitions, making possible to have 
excitations in the forbidden region < uj < 2/z, as ob- 
served experimentally.'^ In this section, we are interested 
in studying the effect on the optical spectrum of doped 
graphene of the different kinds of disorder, considered in 
the previous section. 

In Fig. [5] we compare the numerical results for 
the optical conductivity of doped graphene, consider- 
ing four different types of non-correlated disorder (ran- 
dom potentials, random hoppings, vacancies and hydro- 
gen adatoms) as well as clean graphene. First, one notice 
that the effect of doping is not relevant in the high energy 
part of the spectrum (w 3> /i), and (j(w) follows the same 
behavior discussed in Sees. |III| and [rv] with a peak corre- 
sponding to particle-hole inter-band transitions between 
states of the Van Hove singularities at w ~ 2t. However, 
the spectrum changes dramatically in the infra-red re- 
gion, as shown in the insets of Fig. [5] Therefore, from 
now on we will focus our interest on the effect of disor- 
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FIG. 5. (Color online) Simulation results of the optical con- 
ductivity of doped graphene with diffenrent kinds of non- 
correlated disorders. The chemical potential is fi = O.lt in 
(a) and 0.2f in (b). 



ated to transitions between the above discussed impurity 
band and states at the Fermi level. 

To gain more insight about the effect of disorder in the 
optical conductivity of doped graphene, in Fig. [6]we show 
<j{uj) for different values of /i at fixed concentration of im- 
purities (upper panels), and <7{io) for different concentra- 
tions of impurities and fixed /i. In the first case, the main 
feature is that the conductivity increases as the doping 
decreases, in a qualitative agreement with the experimen- 
tal results.'^ When the chemical potential is fixed and the 
concentration of impurities changes (bottom panels) , one 
observe that the conductivity in the region < a; < 2/i 
grows with na;(j) from cr(w) = for a clean sample to 
<j{uj) « 0.4(To for the larger concentration of impurities 
consider {nx(i) = 0.5%). If we compare to recent ex- 
periments, we notice that a 0.25% of resonant impurities 
would lead to a background contribution similar to the 
one reported by Li et al. for graphene on Si02 ^ whereas 
only a ~ 0.1% of resonant impurities would be neces- 
sary to quantitatively reproduce the results of Chen et al. 
for graphene doped with a high-capacitance ion-gel gate 
dielectric. Finally, we can see that similar results are 
obtained when a sample with correlated on-site potential 
disorder distributed in the form of Gaussian clusters, as 
shown in Fig. [7] Therefore, we conclude that there are 
several kinds of disorder (resonant scatterers and corre- 
lated impurities) that can induce a finite conductivity in 
the infra-red region of the spectrum, as observed exper- 
imentally. It is the whole set of data on DC and AC 
transport from which one may infer the dominant type 
of defects in real graphene. 



der on this low energy part of the spectrum. First, one 
notices that for all kinds of disorder, there is a peak in 
(t(cli) close to a; = 0, whereas at slightly higher energies, 
(t(cli) drops to almost zero for the case of non-resonant 
scatterers (red and green curves), while there is still a 
non-zero background contribution when resonant scat- 
terers are considered (light and dark blue curves). This 
can be understood as follows: for all the cases, disor- 
der leads to a broadening of the bands, which allows for 
intra-band transitions between surrounding states of the 
Fermi level. However, we have seen that resonant impu- 
rities create an impurity band at the Dirac point, with 
the corresponding peak in the DOS at E = 0, whereas 
non-resonant impurities are not so effective in creating 
midgap states. Therefore, the background contribution 
that we find in Fig. [5ja)-(b) between < oj < 2fi for 
samples with resonant scatterers are due to transitions 
between the newly formed impurity band and the con- 
duction band. Taking into account that resonant impu- 
rities are believed to be the main source of scattering in 
graphene,'^^!'^ our results suggest that this kind of im- 
purities could be behind the background contrib ution to 
the optical conductivity observed expcrimcntally.'^'^Eillll 
Finally, notice that the peak observed in (j{uj) for the 
case of resonant impurities at the energy uj ^ fi is associ- 



VI. CONCLUSION AND DISCUSSION 

We have presented a detailed theoretical study of the 
optical conductivity of graphene with different kinds of 
disorder, as resonant impurities, random distribution 
of on-site potential or random renormalization of the 
nearest-neighbor hopping parameter (which can account 
for the effect of substitutional defects). Furthermore, we 
have consider the possibility for the impurities to be cor- 
related or non-correlated. 

For all types of disorder considered, the high energy 
peak at w « 2i, due to inter-band excitations between 
states of the Van Hove singularities of the valence and 
the conduction bands, are always sensitive to disorder, 
getting smeared out proportionally to the strength of dis- 
order. On the other hand, the low energy part of the op- 
tical spectrum is strongly dependent on the type of disor- 
der, as well as its strength and concentration. In general, 
for undoped graphene and in the presence of small dis- 
order of the on-site potentials or in the nearest-neighbor 
hopping between the carbon atoms, the characteristics of 
the single-particle Dirac cone approximation are clearly 
present in the spectrum, and a{ui) « ctq at energies for 
which the continuum approximation applies. This is also 
true when we consider Gaussian hopping parameters. On 
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FIG. 6. (Color online) Optical conductivity of doped graphene with resonant scatterers. Upper panels: Fixed concentration 
of impurities, cr{ijj) for different values of ^. Lower panels: Fixed chemical potential jj,, a{oj) for different concentration of 
impurities. 



the other hand, if there are long-range Gaussian poten- 
tials, the local shifts of the Dirac points leads to electron- 
hole puddles and to the emergence of states in the vic- 
inal region of the Dirac points. As a consequence, we 
observe an enhancement of the optical conductivity in 
the infra-red part of the spectrum. Interestingly, in the 
presence of resonant impurities (vacancies or hydrogen 
adatoms) there appear midgap states which are quasilo- 
calized around the impurities, the number of which is 
proportional to the number of carbon atoms which are 
located around the impurities. Completely random dis- 
tributed (non-correlated) resonant impurities lead to the 
strongest enhancement of zero modes (seen as a promi- 
nent peak in the DOS at zero energy) and also the largest 
effect on the optical spectrum. In fact, for a large enough 
amount of resonant impurities, we obtain a new peak in 
the optical conductivity at an energy w ~ which is as- 
sociated to optical transitions between the midgap states 
and states of the Van Hove singularities. When, for a 
given concentration of impurities, they merge together 
forming clusters, instead of staying uncorrelated, the in- 
fluence of disorder on the electronic properties becomes 
smaller, especially if these clusters form large islands. 



Finally, we have considered the effect of doping on the 
spectrum. Whereas for clean graphene, only inter-band 
processes with an energy larger than uj = 2^ are optically 
active, the presence of disorder leads to a low energy peak 
in a{uj) (associated to transitions near the Fermi level) 
plus a possible spectral weight in the region < ut < 2fi 
for disorders that can create an impurity band at zero 
energy. Most importantly, we have found that a small 
amount of resonant impurities ^ 0.1 — 0.2%, leads to 
a background contribution to cr(a;) between < lo < 
2/x in qualitative and quantitative agreement with recent 
spectroscopy measurements. 
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FIG. 7. (Color online) Same as Fig. [6] but for doped graphene with on-site potential disorder distributed in gaussian clusters. 
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